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We study the exact nuclear time-dependent potential energy surface (TDPES) for laser-induced 
electron localization with a view to eventually developing a mixed quantum-classical dynamics 
method for strong-held processes. The TDPES is defined within the framework of the exact factoriza¬ 
tion [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105,123002 (2010)] and contains the 
exact effect of the couplings to the electronic subsystem and to any external fields within a scalar po¬ 
tential. We compare its features with those of the quasistatic potential energy surfaces (QSPES) often 
used to analyse strong-field processes. We show that the gauge-independent component of the TDPES 
has a mean-field-like character very close to the density-weighted average of the QSPESs. Oscillations 
in this component are smoothened out by the gauge-dependent component, and both components are 
needed to yield the correct force on the nuclei. Once the localization begins to set in, the gradient of 
the exact TDPES tracks one QSPES and then switches to the other, similar to the description provided 
by surface-hopping between QSPESs. We show that evolving an ensemble of classical nuclear tra¬ 
jectories on the exact TDPES accurately reproduces the exact dynamics. This study suggests that the 
mixed quantum-classical dynamics scheme based on evolving multiple classical nuclear trajectories 
on the exact TDPES will be a novel and useful method to simulate strong field processes. 


PACS numbers: 31.15-p, 31.50.-x, 32.80.-t, 33.80.-b, 42.50.Hz, 82.20.-w 


I. INTRODUCTION 

With the advent of attosecond technology ||T]-0, the 
experimentally accessible time-scale has sliifted to that 
of elecfronic motion. If allows fhe observafion of elec- 
fronic motion in real-fime, and even offers fhe confrol 
of election mofion and localizafion via lasers. Several 
groups IZUll have demonsfrated fhaf if is possible fo 
confrol elecfronic mofion in a dissociating molecule and 
localize if selectively on one of fhe producfs of dissoci- 
afion, wifh several differenf sfrafegies. One fechnique 
employs fhe carrier envelope phase (CEP) of a single 
few-cycle laser pulse ITSHSSli , and anofher employs fhe 
fime-delay between fwo coherenf ulfrashorf pulses ||3- 

dzi. 

These experimenfs so far freaf small sysfems (such as 
H 2 and D 2 ), wifh fhe aim of undersfanding the mech¬ 
anisms of localizafion, before applying fhe fechniques 
fo fhe confrol of larger sysfems ETl I28l . Theoretical 
sfudies have a dual role Il7l- l26l : (i) fo help undersfand 
fhe complex correlation between the electron d 5 mam- 
ics and nuclear d 5 mamics, and (ii) to establish methods, 
generally extendable to larger systems, that accurately 
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simulate the coupled electron-nuclear d 5 mamics. Eor 
systems with more than two or three degrees of free¬ 
dom, we must rely on approximate methods, and usu¬ 
ally some kind of mixed quantum-classical approach is 
appropriate, where the electrons are treated quantum- 
mechanically coupled to nuclei described via classical 
trajectories Il29ti32l . Different mixed quantum-classical 
schemes such as Ehrenfest and surface-hopping Il33l - 
l35l , differ in their treatment of the classical nuclear mo¬ 
tion, but use the same form for the potential acting on 
the electrons. Eor d 5 mamics in strong fields, a surface¬ 
hopping scheme between quasi-static potential energy 
surfaces (QSPES) was introduced ll36H38li , and in fact 
applied to the electron-localization problem IITSli . Al¬ 
though this semiclassical approach was shown to repro¬ 
duce the experimental asymmetries reasonably well, it 
is not altogether clear why surface-hopping should give 
good predictions, given its problems associated with 
over-coherence Il39ti4^ . 

In this paper, we will study the possibility of us¬ 
ing a potential derived from first-principles, the time- 
dependent potential energy surface (TDPES) Il47l |48] , 
in a mixed quantum-classical description of the cou¬ 
pled d 3 mamics. This potential arises out of the ex¬ 
act factorization framework where a time-dependent 
Schrodinger equation (TDSE) for the nuclei alone can 
be formulated. The potentials appearing in this equa¬ 
tion capture exactly all coupling to the electronic system 
as well as any external fields, and the resulting nuclear 
wavefunction reproduces the exact nuclear dynamics. 
The scalar potential is denoted the TDPES, and in many 



situations, including all one-dimension problems, the 
TDPES is the only potential acting on the nuclear sub¬ 
system; its gradient therefore yields the exact force on 
the nuclei. For this reason, it is important to gain an un¬ 
derstanding of ifs strucfure, to address both points (i) 
and (ii) above. Therefore, our aim in this paper is to 
find fhe exacf TDPES for the problem of laser-induced 
elecfron localization in a one-dimensional model of Hj, 
compare its structure with potential surfaces more tradi- 
fionally used for strong-field d 5 mamics, and study clas¬ 
sical nuclear dynamics on the exact TDPES with a view 
to developing mixed quantum-classical schemes based 
on the exact factorization. 

Previous work Il49l - I^ has analysed the structure of 
the exact TDPES for a case of field-free d 5 mamics, non- 
adiabafic charge-fransfer in the Shin-Metiu model 1521 , 
finding fhaf much intuition is gained by analysing it 
in term of fhe Bom-Oppenheimer (BO) potenfial en¬ 
ergy surfaces (BOPESs), and fhaf such an analysis en¬ 
ables connections fo be made wifh fradifional approxi¬ 
mate methods for coupled elecfron-ion d 5 mamics, such 
as surface-hopping. Further, it was found that evolving 
an ensemble of classical nuclear trajectories on the exact 
TDPES accurately reproduces the exact nuclear d 5 mam- 
ics ISTl . 

We will show here that analogous conclusions can be 
drawn for the laser-induced electron localization prob¬ 
lem: an ensemble of classical nuclear trajectories evolv¬ 
ing on the exact TDPES accurately reproduces the exact 
nuclear d 5 mamics, and analysis in terms of the QSPESs, 
which play the role of the BOPESs when strong fields 
are presenf, is helpful. The TDPES nafurally separates 
info a gauge-independenf parf and a gauge-dependenf 
parf. We show thaf the density-weighted average of 
fhe QSPESs approximates the gauge-independent com¬ 
ponent, which is rather oscillatory and the force on the 
nuclei resulting from ifs gradient is incorrect. Once the 
gauge-dependent component of fhe TDPES is included, 
fhe oscillations smoothen ouf: togefher, fhey yield the 
correct force on fhe nuclei. Further, we find fhat, once 
localization begins to set in, the gradient of the exact TD¬ 
PES at the location of the mean nuclear position, tracks 
that of one QSPES and fhen swifches to the other, resem¬ 
bling the picture provided by the semiclassical surface¬ 
hopping approach IIT3l 15711551 . 

A multiple trajectory Ehrenfest d 5 mamics simulation 
shows that although the nuclear d 5 mamics is reasonably 
reproduced, an incorrect electron localization asymme¬ 
try is obtained. The error can be related to the incor¬ 
rect BO projections of the electronic wavefunction. The 
fact that the Ehrenfest d 3 mamics yields inaccurate elec¬ 
tron d 5 mamics can be anticipated from our recenf work 
on the exact electronic-TDFES 1551 : in this complemen¬ 
tary picture, instead of asking what is the exact potential 
acting on nuclei in an exact TDSE for nuclei, one asks 
what is the exact potential acting on electrons in an ex¬ 
act TDSE for fhe electronic subsystem. We found 1551 
fhaf the exact electronic-TDPES is significantly differ¬ 


ent from fhe pofenfial acting on electrons in the usual 
mixed quantum-classical schemes - including Ehrenfest 
as well as surface-hopping schemes - yielding signifi¬ 
cant errors in the prediction of fhe elecfron localizafion 
asymmefry. The results of fhe present paper suggest 
that, instead, mixed quantum-classical schemes based 
on evolving multiple classical trajectories on the exact 
TDPES (or good approximations to it) will be a useful 
mefhod fo simulate strong field processes. 

This paper is organized as follows. In section II, we 
review two different concepts of pofential energy sur¬ 
faces for TD processes in laser fields: the QSPES and the 
exact TDPES. In section III we compare the features of 
these potentials for electron localization d 5 mamics in the 
dissociation of a model molecule induced by time- 
delayed coherent ultra shortlaser pulses. We show the 
exact TDPES gives the correct force acting on nuclei, so 
evolving multiple classical trajectories on it reproduces 
the correct nuclear wavepacket d 5 mamics. The force ob- 
fained from surface-hopping between QSPESs can ap¬ 
proximately reproduce such an exact force once local¬ 
ization begins to set in. We also compute multiple tra¬ 
jectory Ehrenfest dynamics and reveal how it fails to re¬ 
produce electron localization d 5 mamics while it reason¬ 
ably reproduces the nuclear d 3 mamics. In section IV we 
summarize the results and remark on the future direc¬ 
tions. 


II. THEORY 

A. Quasi-static potential energy surface 

In this section we first review the concept of the 
QSPES introduced for fhe description of molecules in 
strong-fields. The QSPES has been thoroughly dis¬ 
cussed in earlier works IIT511561158115411601 , but we here 
give a discussion particularly relevant for the electron 
localization d 5 mamics problem in the dissociation of 
H+. 

For fhis problem, the essential physics is contained 
in the two lowest field-free electronic states of fhe BQ 
Hamiltonian, i.e., the Iscr^ and 2p(T„ states, and the full 
molecular wavefunction r, t) of the system can be 
expressed as 

T'(i?, r, t) = Xg{R, + Xu{R, . (1) 

Here Xg(f?,f) and XuiR,t) describe nuclear wavefunc- 
tions that exist in the IscTg and 2ptT„ states respectively, 
functions of fhe intemuclear disfance R and time t, 
and ^fj(r) and (/)^(r) describe fhe Iscg and 2pcr„ elec¬ 
tronic wavefunction respectively, which parametrically 
depend on R. Since (()fj(r) and (/>^(r) are bonding and 
anfi-bonding combination of Is atomic orbitals, a coher¬ 
ent superposition of them provides the localized elec¬ 
tronic states = ;^(</'fi(r) ± that have 

the electron on either the left or the right proton. These 
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states form a convenient basis in which to monitor the 
electron localization asymmetry. In the experiment, in¬ 
teractions of fhe molecule wifh fhe fime-delayed infra¬ 
red laser field in fhe course of fhe dissociafion provides 
a coupling of </>fj(r) and creafing a coherenf su- 

perposifion sfafe, and, insfead of Eq. if is insfrucfive 
fo wrife: 

( 2 ) 

where xieft(i?, t) and Xright(f?, f) are defined as fhe nu¬ 
clear wavefuncfions fhaf exisf in connecfion wifh (r) 
and Measuremenfs of ion fragmenf asymme- 

fries leff or righf along fhe polarizafion axis direcfly re- 
lafe fo x\eft{R, t) and Xright(i?, t). 

While fhe field-free sfafes above are useful fo analyse 
fhe asymmefry, fo undersfand fhe fime-developmenf of 
fhe localizafion if is helpful fo consider a fhird, fime- 

dependenf, basis, fhe TD quasisfafic sfafes, f), 

also known as phase-adiabafic sfafes. These sfafes are 
defined as insfanfaneous eigensfafes of fhe insfanfa- 
neous elecfronic Hamilfonian t), defined by 

-Ff}laser(r,f), (3) 

i.e., 

= S^^\R,t)cl>f^^{v,t) (4) 

where (i?, f) are fhe quasisfafic pofenfial energy 
surfaces (QSPESs). Wifhin our fwo-sfafe model we may 
wrife 

= cW(i?,f)0fj(r) -hcW(i?,f)0^(r), (5) 

so fhaf fhe eQ®b) of Eq. § are given by fhe eigen¬ 
value equafion: 

( \ ^^QS(i) ( cf 

\ ) \ ) \ c« 

( 6 ) 

Therefore we can express fhe QSPESs in ferms of fhe 
BOPESs e®oW(i?) as 


The molecular wavefuncfion expressed in ferms of 
quasi-sfafic sfafes is 

T'(i?,r,f) = xT + xT 

( 10 ) 

Nofe fhaf fhe nuclear wavefuncfions (i?, t) and 
X 2 ^(i?, f) thaf are connecfed fo fhe quasi-sfafic sfafes 
j^QS(i) (r, t) can be expressed in ferms of 

Xleft f) and Xright 

^ [xieft (f?, t) (cos 9 +sin 9) 

+ Xright (cos 6»-sin 6»)] 

[xieft iR,t){-cos9 + sin 9) 

+ Xright (R, t) (cos 6> -h sin 6»)]. 

which can be used fo exfracf fhe elecfron localizafion 
from x?^(f?, t) and x^^(f?, t). 

A semi-classical surface-hopping model based on 
QSPESs has recenfly been ufilized fo undersfand and 
reproduce fhe elecfron localizafion dynamics and asym¬ 
mefry imisziiMi in H J. In fhis approach, an ensemble 
of classical nuclear frajecfories evolve on one QSPES or 
fhe ofher QSPES, making insfanfaneous hops befween 
fhem as defermined by a Landau-Zener formula. If 
was shown fhaf fhe elecfron localizafion sefs in a region 
where fhe d 5 mamics is infermediafe befween adiabafic 
and diabafic: fhe ensemble of nuclear frajecfories fra- 
verses several laser-induced avoided crossings befween 
fhe QSPESs. This semi-classical mefhod gives asym¬ 
mefry paramefers in reasonably good overall agreemenf 
wifh fhaf obfained from fhe full TDSE alfhough fhe de¬ 
fails differ. The agreemenf lends some hope fo fhe use of 
fhis semiclassical scheme fo simulafe coupled elecfron- 
ion d 5 mamics in confrol problems in more complicafed 
sysfems; however, af fhe same fime a furfher under- 
sfanding of fhe errors in fhe defails is desirable. We will 
analyse fhis approach by comparing fhe QSPESs wifh 
fhe exacf TDPES, which we will review in fhe nexf sec- 
fion. 


xf(i?,i) = 
X2^{R,t) = 


gQ®(1.2) ^£B0(1.2) gBO(2,l) gj^2 

± ('('fjlwiaserl'/'^) Sm29{R,t) 


B. Exact time-dependent potential energy surface 


(7) 

and the electronic quasi-static eigenstates in terms of fhe 
BQ sfafes, 

= cos9{R,t)(l)%ir) + sm9{R,t)(j)'^{r) 

= sin6»(i?,f)(/)fj(r) - cos9{R,t)(j)'^{r), 

where fhe TD mixing paramefer 9{R, t) is given by 


In Ref. Il47l |48l , if was shown fhaf fhe full molecular 
wavefuncfion 4^(r, R, t) which solves fhe TDSE 


fT'I'(r, R, f) = zd(T'(r, R, f) (12) 

can be exacfly factorized fo fhe single producf 

(i, S, t) = x(E, f )%(i, t) (13) 


tan 29{R, t) 


eBO(l)(^)_gBO(2)(^)' 


of fhe nuclear wavefuncfion x(B;jf) and fhe elecfronic 
wavefuncfion $R(r, f) fhaf paramefrically depends on 
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the nuclear positions R and satisfies the partial normal¬ 
ization condition 


dr|$R(r,t)p = 1 VR,h 


(14) 


Here, the complete molecular Hamiltonian is 

H = f^m + (15) 

and Hbo (Hj R) is the BO electronic Hamiltonian, 

HBO=fe (l) + Ihee (g) + Wen (l, E) + Wnn{R)- (16) 

Note that fn = - 2 l& and R = - ^ 

are the nuclear and electronic kinetic energy operators. 
Wee, Wen and Wnn are the electron-electron, electron- 
nuclear and nuclear-nuclear interaction, and (R, t) 
and f)|^((r, t) are time-dependent (TD) external poten¬ 
tials acting on the nuclei and electrons, respectively. 
Throughout this paper R and r collectively represent 
the nuclear and electronic coordinates respectively and 
h=i. 

Returning to Eq. | (T^ , the stationary variations of fhe 
quanfum mechanicaTacfion wifh respecf fo $R(r, t) and 


E ^ - o.. + t) + e(R, t) 


,a—l 


2M„ 


x(R,t) 


= idtx{R,t)- 

(18) 


Here, e(R, t) is fhe exacf nuclear TDPES 


e(R,f) = (%(f) HBo+vt,t{x,t) + U^e7^-tdt %(f))^, 

(19) 

[/g°“P[<i>R, x] is fhe "elecfron-nuclear coupling opera¬ 
tor", ^ 






a—1 

-*VaX 


[-fVa - AQ(R,f)]' 


( 20 ) 


X 


+ Aq(R, f)j ( —fV C ~ Aa(R, t)) 


and Aq, (R, f) is fhe TD vecfor pofenfial pofenfial. 


Aq (R, f) — (4>R(f) — iVa ‘l’R(f) 


( 21 ) 


The symbol ( • ) j. indicafes an integration over elecfronic 

coordinafes only Note fhaf fhe PNC makes fhe fac¬ 
torization iTTsJ unique up to within a (R, f)-dependent 
gauge transformation. 


x{R,t) ^ x{R,t) = 




x{R,t) 




( 22 ) 


x(R, f) under the condition (14i lead to the following 
equafions of mofion for x(Ej ijand 4)R(r, t): 

{HboR. R) + i;L,(p, t) + X] - eiR, t)) $g(£, t) 

(17) 


and Eqs. (TT} and ( [TS) are form invarianf under fhis 
fransformafion while the scalar potential and the vector 
potential transform as 

e(R,f) = e(R,t)+dt0(R,f) (23) 

Aa{R,t) = Aa{R,t) +VaO{R,t). (24) 

The equation for the exact nuclear wavefunction, 
Eq. (18 1 , is Schrodinger-like, and the TD vector poten¬ 
tial (211 and TD scalar potential ( [T^ that appear in it, 
exactly govern the nuclear d 5 mamics. It is important 
to note that x(Rj t) can be interpreted as the exact nu¬ 
clear wave-function since it leads to an A-body nu¬ 
clear density, r(R, t) = |x(E, f)P, and an A-body cur¬ 
rent density, Jq(R, t) = Im(x*(R, f)Vax(R,f)) + 

r(R, f)AQ,(R, t) , which reproduce the true nuclear A- 

body density and current density Il48l obtained from fhe 
full wave-function 'I'(r, R, t). 

In our previous work fhe shape of fhis exacf TDPES 
has been useful fo inferpref d 5 mamics for bofh a sfrong 
field process (sfrong-field dissociation of H^) Il47l l48l as 
well as for field-free d 5 mamics of non-adiabatic charge- 
fransfer Il49ll^ . In particular, in the field-free case, a 
defailed sfudy of fhe form of ifs gauge-dependenf and 
gauge-independenf parfs proved insfrucfive fo under- 
sfand ifs effecf on fhe nuclear d 3 mamics, and fhe sfruc- 
fure fo be expecfed for general field-free problems. Im- 
porfanfly, in a mixed quanfum-classical description, the 
gradient of fhis exacf TDPES gives uniquely fhe correcf 
force on fhe nuclei, and if was shown, in the field-free 
case, fhaf an ensemble of classical frajecfories evolving 
on fhe exacf TDPES accurately reproduces fhe exacf nu¬ 
clear wavepackef d 3 mamics. We now consider a defailed 
sfudy of fhe form of fhe exacf TDPES for fhe presenf case 
of d 5 mamics in external fields, wifh fhe aims of address¬ 
ing fhree quesfions. Eirsf, does running classical nuclear 
d 5 mamics on fhe exacf TDPES reproduce fhe d 5 mamics 
of laser-induced elecfron localization? Second, how are 
fhe QSPESs related to the exact TDPES? Third, can we 
see hints of fhe semiclassical surface-hopping mefhod 
in fhe exacf TDPES? 


III. RESULTS AND DISCUSSION 
A. Theoretical model 

We employ a one-dimensional model of fhe H^ 
molecule fo sfudy elecfron localization d 3 mamics 
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achieved by time-delayed coherent ultra short laser 
pulses IZl iSl [T3l . In the experiment, first an ultravio¬ 
let (UV) pulse excites to the dissociative 2pcrt, state 
while a second time-delayed infrared (IR) pulse induces 
elecfron fransfer befween fhe dissociafing atoms. In our 
model, we sfarf fhe d 5 mamics after fhe excifation by fhe 
UV pulse: fhe wavepackef sfarfs af f = 0 on fhe firsf 
excifed sfafe (2p(T„ sfafe) of as a Frank-Condon pro- 
jecfion of fhe wavefuncfion of fhe ground sfafe, and fhen 
is exposed fo fhe IR laser pulse. The full Hamiltonian of 
fhe system is given by 

H{R,z,t)=f^{R) + W^\z,t) 

= Tn{R) + Te{z) + Wnn{R) + Wen{z, R) + '0laser(2:, t) 

(25) 


where R is fhe infemuclear disfance and z is fhe elec- 
fronic coordinate as measured from fhe nuclear cen- 
fer of mass. The kinetic energy ferms are T„(i?) = 

“ 2)b ^ and, fe{z) = - ^ respectively, where fhe 
reduced mass of fhe nuclei is given by /i„ = 
and reduced elecfronic mass is given by pe = 2 Mh+i 
(Mh is fhe profon mass). The inferacfions are soff- 


Coulomb: 

1 


W„n{R) = 


\/0.03+R^ 


, and Weniz,R) = 

- / „ = — , ^5= (and Wee = 0). The 

yi.o+(z-f)2 ^i.o-r(2+f)2 ^ 

TR pulse is described wifhin fhe dipole approxima¬ 
tion and lengfh gauge, as vl^^{z,t) — E{t)qeZ, where 

E{t) = Ugexp — cos(a;(f — At)), and fhe re¬ 


duced charge Qe = 2 Mh+i ' wavelengfh is 800 nm 
and fhe peak intensify Iq = Eq = 3.0 x lO^^W/cm^. The 
pulse duration is r = 4.8/s and At is fhe fime delay be¬ 
tween the UV and IR pulses. Here we show the results 
of Af = 7 fs. 

We propagate the full TDSE 


H{z, R,t)'^{z, R,t) = R,t) (26) 


numerically exactly to obtain the full molecu¬ 
lar wavefuncfion 'i>{z,R,t), and from it we cal¬ 
culate the probabilities of directional localiza¬ 
tion of the electron, P±, which are defined as 
P+{-) = /2>(<)o^)l^' These are shown 
as the green solid (P-) and red dashed (P+) lines in 
Fig. [^. It is evident from this figure that considerable 
electron localization occurs, with the electron density 
predominantly localized on the left (negative z-axis). 

Furthermore, we calculate the population d 3 mamics 
of the BO states (t>%{z) (green solid) and (j^Ri^) (red 
dashed) (Fig.j^) during dissociation, as well as the pop¬ 
ulation d 5 mamics on the 1st quasi-static state {z, t) 
(green solid) and 2nd quasi-static state (j}^^^'^\z,t) (red 
dashed) (Fig. m the relative simplicity of the lat¬ 
ter demonstrate the usefulness of the QS basis for 
laser-induced processes. We then plot the QSPESs 
e^^^^'>{{R{t)),t) (green solid) and e^^^‘^\{R{t)),t) (red 
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FIG. 1. (a) 4.8 fs FWHM 800 nm laser pulse, (b) Elec¬ 

tron localization probabilities along the negative (green solid 
line) and the positive z-axis (red dashed line) as a function 
of time, (c) Population dynamics during dissociation on the 
BO state 4>%{z) (green solid) and (pniz) (red dashed), (d) 
Population dynamics during dissociation on the 1st quasi¬ 
static state (green solid) and 2nd quasi-static state 

^QS( 2 ) dashed), (e) Quasi-static potential energy sur¬ 

faces e^iSCi) (green solid) and (R, t) (red dashed) for 
a nuclear trajectory {R) (t) that tracks the expectation value of 
the internudear distance. The blue curve shows the transi¬ 
tion probability given by a Landau-Zener formula (Eq. 18 of 
Ref. 


dashed) evaluated at a nuclear trajectory {R{t)) that 
tracks the expectation value of the internudear distance. 
These results coincide qualitatively with the previous re¬ 
sults reported by Kelkensberg et al. IIT^ Panels b, d, and 
e, suggest that the electron localization is determined 
by the passage of the dissociating molecule through a 
regime where the laser-molecule interaction is neither 
diabatic nor adiabatic. As discussed in the previous sec¬ 
tion, the semiclassical scheme, with the avoided cross¬ 
ings between the QSPES inducing the trajectories to hop 
between them, reproduces the general behavior. Next, 
we will compare the exact TDPES with the QSPES to un¬ 
derstand the relation between the two, shed some light 
on the surface-hopping scheme, and find the exact force 
on classical nuclei. 
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B. Exact TDPES vs. QSPES 


First we show the exact TDPES for this process in 
Fig-H We calculate the TDPES in the gauge where 
the vector potential A{R,t) is zero l48l , so the TD¬ 
PES €{R,t) is the only potential acting on the nuclear 
subsystem. It is instructive to express the TDPES as 
the sum of fhe gauge-independenf ferm egi{R,t) and 
fhe gauge-dependent term €gd{R, t) as done in previous 
studies iSSlIlgl : 

e{R,t) = ^gi{RA) + ^gd{R:t) (27) 

where 

eg,{R,t) = {‘^R{t)\HBO + Glaser + (28) 

and 

egd{R,t) = - idt l^’fl(f))^ ■ (29) 

In Eig. 0 e(i?, f)(black solid), t)(blue solid) and 

egd{R, f)(orange solid) are plotteci at nine different times, 
along with the two lowest BOPESs, e®°(^i(i?) and 
gBO( 2 )(^) (Note that the TDPES e(i?, f)(black solid) 
and its GD component egd{R, f)(orange solid) have been 
rigidly shifted along the energy axis). 

We also plot the exact nuclear density |x(i?,f)p (green 
solid) and the nuclear density reconstructed from evolv¬ 
ing an ensemble of 800 classical frajecfories on fhe ex- 
acf TDPES (red dashed) ISTI af each time. The closeness 
of fhese lasf fwo curves shows fhaf a mixed quanfum- 
classical scheme for the electron localization process is 
appropriate and that the exact TDPES e{R,t) gives the 
correct force acting on classical nuclei in such a scheme. 

In previous work 1149115X1 , sfep-like feafures of Cgi {R, t) 
and egd{R,t) in fhe field-free non-adiabafic process in 
fhe vicinify of fhe avoided crossing have been shown. 
In particular, after passage through the avoided cross¬ 
ing, where the nuclear wavepacket had spatially sepa¬ 
rated on two BOPESs, the GI component tracked one BO 
surface or fhe ofher, wifh a sfep befween fhem, while 
fhe GD componenf was piecewise flaf, buf wifh a sfep 
in fhe same region wifh opposife sign. The nef TD¬ 
PES was overall more smoofh fhan eifher of fhe compo- 
nenfs. Here, we find again very inferesfing feafures of 
egi{R, t) and €gdiR, t). Eirsf nofe fhaf bofh egi{R, t) and 
egd{R,t) shows many small hills and valleys affer fhe 
laser-induced nonadiabatic transifions begin, buf wifh 
opposife slopes fo each ofher, so fhaf fhese sfrucfures 
largely cancel each ofher when fhe exacf TDPES e{R, t) is 
consfrucfed (much like fhe near-cancellafion of fhe sfeps 
in fhe field-free case). Like fhe field-free case, bofh fhe 
GI and GD ferms are imporfant fo consider fo predicf fhe 
correcf nuclear dynamics. Second, in fhe presenf sfrong- 
field case, unlike fhe field-free examples sfudied in 1491 - 
ISTl , egi{R, t) does nof piecewise frack one BOPES or fhe 
ofher. However, if does frack a density-weighted QSPES, 
as we will show nexf. 



FIG. 2. Snapshot of the exact TDPES e{R, t) (black solid) 
, its gauge-invariant part tgi{R,t) (blue solid) and gauge- 
dependent part egd{R,t) (orange solid) at indicated times 
along with two lowest BOPESs (black dashed). Furthermore, 
the exact nuclear density \x{R,t)\'^ (green solid) and the nu¬ 
clear density reconstructed from the multiple trajectory dy¬ 
namics on the exact TDPES (red dashed) for each time are also 
plotted. 


In Eig.l^ we show e(i?, f)(black solid) (which is again 
rigidly shifted along the energy axis) and the gauge- 
invariant part egi{R,t) (blue solid) together with the 
QSPESs e^^^R{R,t) (green solid) and e^^^‘^'>{R,t) (red 
solid). We find fhaf fhe oscillations in fhe gauge- 
invariant part of exacf TDPES Cgi {R, t) (blue solid) 
fend to step between the two QSPESs: |x^®(i?,f)p 
and |x^^(f?, f)P are also plotted in Eig. and we see 
that egi{R,t) tends towards the QSPES whose popu¬ 
lation is dominant, i.e. when |x5^®(Ff, t)P is larger 
than |x?^(f?i f-gi{R, t) approaches to t) and 

when |x?^(^if)P is larger than |x?^(f?, t)P egi{R,t) ap¬ 
proaches to e^®i^^(i?,f). In fact, egi{R,t) lies practically 
on top of fhe fhe weighfed average of fhe quasi-sfatic 
surfaces e2fe(^> 0- 


eTeiR^t) 


\xTiR,tr 

|xf(i?,f)P + |xf(i?,f)P 

\xT(R,tr 

|x?"(i?,f)P + |xf(i?,f)P 




(30) 


This is ploffed wifh lighf blue line in Eig. Therefore 
fhe weighfed-average of fhe QSPESs approximafes fhe 
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FIG. 3. Snapshots of the gauge-invariant part of exact TO¬ 
PES egi{R,t) (blue solid), QSPESs f) (green solid) 

and (red solid), and the weighted average of 

the QSPESs e^^^{R,t) (light blue solid) at indicated times. 
|x?^(f?, f)P (green) and 1)1^ (red) are also plotted. 


gauge-invariant part of exact TDPES egi{R,t), but not 
the full exact TDPES e{R,t). In fact, this is quite anal¬ 
ogous to the previous results on the field-free passage 
through an avoided crossing ||49H5T]| : there, at the times 
considered, the density-weighted average collapsed to 
one BO surface or the other except in the intermedi¬ 
ate (step) region, because the spatial separation of the 
parts of the density projected onto different BO surfaces 
meant that in the field-free analog to Eq. 1^ the prefac¬ 
tors of each term was either one or zero. Flere it is evi¬ 
dent that the density does not spatially separate (Eig. |^, 
i.e. the projections on to the QSPESs overlap. One can 
make entirely analogous statements in both cases: the 
density-weighted average of the BOPES approximates 
the gauge-invariant part of exact TDPES egi{R, t) in the 
field-free case, and the density-weighted average of the 
QSPES approximates the gauge-invariant part of exact 


TDPES egi{R, t) in the presence of strong fields. 

To confirm the relationship between egi{R,t) and 
t), we consider the expansion of the complete 
wavefunction with the two lowest quasi-static states 
(Eq. [To) . Then the exact electronic conditional wave- 
function ^r{z, t) is expressed as: 








{z, t) -\- 


xiR,t) 


aQS(2) 

Vr 


{z,t). 


(31) 


Then we realize: 

_ lx? QS(1) I lx? QS(2) ( 32 ) 

\x{R,tW \x{R,tW 

= eTeiR,t). 


Since egi{R,t) = {^r{zR)\H^‘^ + v\asev\'iR{z,t))^ + 
2 ^($i^(z,^)|(-^^ - A{R,t)f\^R{zR))^, we can con¬ 
clude 


eg^iRR) ^ e‘^iiR,t), (33) 

because 0{M~^) term gives a much smaller contribu¬ 
tion. 

To reproduce the correct d 5 mamics, however the ef¬ 
fect of egd{R, t) is crucial to include, as in the field-free 
case studied before lai. In the gauge we have cho¬ 
sen A{R^ t) = 0, but we note that if instead we choose 
the gauge where CgdiR, t) = 0 then the vector potential 
A{R, t) will be non-zero, and will be responsible for the 
role of effectively reducing the oscillatory structure in 
the GI term. 

In Eig. [^ we plot the gradient of the different poten¬ 
tials computed on the trajectory of mean nuclear dis¬ 
tance {R){t), as a more direct probe of the force on the 
nuclei. The black line, which is the gradient of the ex¬ 
act TDPES ^i{{R{t))), gives the exact force on the nu¬ 
clei. Eirst we immediately notice that the gradient of 
the weighted average of the two QSPES ^eSfe((l?(0)) 
(light blue line) (equivalently, the GI component (blue 
line)) is completely different from the exact force. A 
semi-classical simulation on the weighted average of the 
two QSPES would not give the correct nuclear d 5 mam- 
ics. We observe instead that, as the localization sets in, 
the exact force ^e((^(f))) coincides with the gradient 
of one or the other QSPES (red or green). This sup¬ 
ports the idea of semiclassical surface-hopping between 
QSPES IIT3l l37l |38] at least after the localization begins 
to set in (time ^6 fs): the exact force on the nuclei is 
given by the gradient of fhe exacf TDPES, and, when 
evaluafed af fhe mean nuclear position, coincides wifh 
the force from one QSPES or fhe ofher, making fransi- 
fions between fhem af fheir avoided crossings. This ex¬ 
plains why the semiclassical simulations of Ref. 1131 had 
a reasonable agreement with the exact results. Eurther- 
more the figure shows fhe imporfanf role of fhe gauge- 
dependenf parf egd{R,t)f wifhouf fhis ferm, fhe force 
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FIG. 4. Time evolution of the gradient of each TDPES at 
position Green line: Red line: 

Blue line: Light blue line: 

Black line: ^e({i?(t))). 


on the nuclei would be more oscillatory and quite dif¬ 
ferent (blue line in the figure). We nofe fhaf if insfead 
we choose fhe gauge where CgdiR, t) — 0 fhen fhe vec¬ 
tor pofenfial A{R,t) will be responsible for fhe role of 
effectively reducing fhe oscillatory sfrucfure in fhe GI 
term. As sfafed above, when we choose fhe gauge where 
egd{R,t) = 0, fhen fhe vector pofenfial A{R,t) plays 
fhe role of if according to fheir relafionship: A{R, t) = 
f* dt' {-dRCgdiR, t' j) ISll. 


C. Multiple trajectory Ehrenfest dynamics 

Given that there are several avoided crossings dur¬ 
ing the localization d 5 mamics, one might ask how well 
a mean-field surface to propagafe fhe elecfrons would 
work. To fhis end, we run a mulfiple-frajecfory Ehren- 
fesf calculation ICTI , and compare fhe elecfron and nu¬ 
clear densities wifh fhe exacf ones. 

In fhe upper panel of Fig. [5 we plof fhe condi¬ 
tional elecfron densify obfained from fhe 

exacf calculafion af indicated times. In fhe lower 
panel, we plof ifs squared expansion coefficienfs of 
fhe Bom-Oppenheimer expansion \Cg{R,t)\^ (green) 
and |G„(i?,f)p (red) (T>i^(z,^) = Cg{R,t)^R{z) + 

G„(i?, ^$^(z)), along wifh fhe nuclear densify (black). 
In Fig.l^ we plotted elecfron densify \^{z,t\Rci{t))\'^ ob¬ 
fained from fhe multiple frajecfory Ehrenfesf d 3 mamics 
calculafion af fhe indicafed times (plotted for all 800 fra- 
jecfories Rci{t)). The lower panel shows fhe squared 
expansion coefficienfs of fhe Born-Oppenheimer ex¬ 
pansion \Cg{Rci{t))\‘^ (green) and |G„(i?ci(f))P (red) 
of fhe elecfronic wave function ^{z,t\Rci{t)) obfained 
from mulfiple frajecfory Ehrenfesf d 3 mamics calculafion 
{^iz,t\Rcl{t)) = Gg(i?ci(t))$fj(z) + G„(i?ei(f))nW)- 



t=12.0fs 



8 10 12 14 

R (a.u.) 


FIG. 5. Upper panel: Conditional electron density \^R{z,t)\’^ 
obtained from the exact calculation at the indicated times. 
Lower panel: Squared expansion coefficients of the Born- 
Oppenheimer expansion \Cg{R,t)\'^ (green) and |G„(i?,t)p 
(red) of the exact conditional electronic wave function ^r{z, t) 
(^R{z,t) = Gg(ii, t)<l?fj(z) + Cu{R,t)^R{z)) at the indicated 
times. The exact nuclear density is also plotted (black). 


We also show the nucler densities reconstructed 
from the distribution of classical trajectories obtained 
from multiple trajectory Ehrenfest d 5 mamics calculation 
(black circle line). 

Comparing the top panels of these figures shows that 
Ehrenfest dynamics gets the overall structure of the elec¬ 
tronic conditional probability reasonably well, however 
not exactly; in fact, these differences lead to an incorrect 
prediction of the localization asymmetry. For example, 
at f = 12.0 fs at the internuclear separation where the 
nuclear density is peaked, the projections shown in the 
lower panels, onto the BO surfaces predicted by Ehren¬ 
fest are each close to 0.5, while the exact are closer to 
0.6 and 0.4. Given the nature of the BO g and u states 
in terms of the left and right basis (Sec. |II A|, this sug¬ 
gests the localization asymmetry predicted by Ehrenfest 
is close to 1:0 while the exact is close to 0.8:0.2. Indeed 
this is verified by the calculation of the asymmetry. Fur¬ 
ther, throughout the width of the nuclear wavepacket, 
the Ehrenfest projections remain close to 0.5, while the 
exact projections fall away, indicating there is a larger 
degree of decoherence in the exact dynamics, missed 
in the Ehrenfest d 5 mamics. The differences in the con¬ 
ditional wavefunction and the BO projections is even 
greater where the nuclear density is small (i? = 8 ^ 10 
and R= 12 14). 

In the field-free problem of non-adiabatic charge- 
transfer MMM, multiple-trajectory Ehrenfest dy¬ 
namics failed, and this might have been expected given 
that the density spatially separates (branches) onto two 
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FIG. 6. Upper panel: Electron density \^{z,t\Rci{t))\’^ 
obtained from multiple trajectory Ehrenfest dynamics cal¬ 
culation at the indicated times (plotted for all trajecto¬ 
ries Rci{t))). Lower panel: Squared expansion coef¬ 
ficients of the Born-Oppenheimer expansion \Cg{Rci{t))\'^ 
(green) and |C'„(i?ci(^))r (red) of fhe elecfronic wave func¬ 
tion ^{z,t\Rci{t)) obtained from multiple trajectory Ehren¬ 
fest dynamics calculation ($(«, t\Rci{t)) = Cg{Rci{t))$^i^{z) + 
Cu{Rci{t))^R{z)) at the indiacted times. Nuclear density re¬ 
constructed from the distribution of classical trajectories are 
also plotted (black circle line). 


different BO surfaces. In the present case, the nuclear 
density does not split in space, and actually predicts the 
nuclear dynamics quite well, but the errors in the elec¬ 
tronic d 5 mamics are more significant. Further, it is the 
same potential that evolves the electrons in the Ehren¬ 
fest calculation as in surface-hopping calculations, and 
this same potential was shown to lack significant struc¬ 
tures that the exact potential acting on the electron sub¬ 
system (e-TDPES) in Ref. 1153 1 has. 


IV. CONCLUSIONS AND OUTLOOK 

The TDPES and vector potential arising from the ex¬ 
act factorization of the molecular wavefunction exactly 


accounts for the coupling to the electronic subsystem 
as well as coupling to external fields and so it is im¬ 
portant to understand their structure, and to relate this 
to the QSPES which is traditionally used, in order to 
be able to develop accurate practical mixed quantum- 
classical methods for strong-field d 5 mamics. In this pa¬ 
per, we have studied the topical phenomenon of laser- 
induced electron localization in the dissociation of Hj, 
choosing a gauge where the TDPES is only potential act¬ 
ing on the nuclear system. We found that the gauge- 
independent component of the TDPES has a mean-field- 
like character very close to the density-weighted aver¬ 
age of the QSPESs and yields an oscillatory force on the 
nuclei. The gauge-dependent component of the TDPES 
smoothens the oscillations of the gauge-independent 
component and together they lead to the correct force. 

We demonstrated that rurming an ensemble of classi¬ 
cal nuclear trajectories on this exact TDPES accurately 
reproduces the exact nuclear d 5 mamics. We found that 
the force obtained by considering surface-hopping tran¬ 
sitions between QSPESs at the laser-induced avoided 
crossing approximates this exact force, after the localiza¬ 
tion begins to set in. We showed that errors in multiple- 
trajectory Ehrenfest d 5 mamics are less significant for the 
nuclear d 3 mamics than for the electronic d 3 mamics ex¬ 
plored in Ref. Il53ll , where it was shown that Ehrenfest 
yields an incorrect electron localization asymmetry. It is 
worth noting that the potential acting on the electrons 
in Ehrenfest dynamics and in surface-hopping schemes 
lack important step and peak features that the exact po¬ 
tential acting on the electronic system (the e-TDPES) has. 
Therefore the results of this study show that to repro¬ 
duce the laser-induced electron localization d 3 mamics 
accurately by means of a mixed quantum-classical dy¬ 
namics scheme, we have to go beyond the traditional 
methods such as surface-hopping or Ehrenfest meth¬ 
ods. Our results here encourage the development of 
mixed quantum-classical schemes based on Eqs ( [T7| and 
| (T^ 1621 to simulate strong-field processes. 
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